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Abstract 

Riemami ellipsoids model rotating galaxies when the galactic velocity field is a linear 
function of the (-artesian coordinates of the galactic masses. In nuclear physics, the kinetic 
energy in the linear velocity field approximation is known as the collective kinetic energy. 
But, the linear approximation neglects intrinsic degrees of freedom associated with nonlinear 
velocity fields. To remove this limitation, the theory of symplectic. dynamical symmetry 
is developed for classical systems. A classical phase space for a self-gravitating symplectic 
system is a co-adjoint orbit ol the noucompact group Sp(3,R). The degenerate co- adjoint 
orbit is the 12 dimensional homogeneous space Sp(3,R)/U(3), where the maximal compact 
subgroup U ( 3 ) is the symmetry group of the harmonic oscillator. The Hamiltonian equations 
of motion on each orbit form a Lax system A = [A r , F], where X and F are elements of the 
symplectic Lie algebra. The elements of the matrix X are the generators of the symplectic 
Lie algebra, viz., the one-body collective quadratic functions of the positions and momenta 
of the galactic masses. The matrix F is composed from the self-gravitating potential energy, 
the angular velocity, and the* hydostatic pressure. Solutions to the Hamiltonian dynamical 
system on Sp(3,H)/L(3) are given by symplectic isospectral deformations. The Casimirs of 
Sp(3,R), equal to the traces ol powers of A . are conserved quantities. 


1 Riemann Ellipsoids 

A remarkably unified pict ure of rotating systems is attained by adopting an algebraic perspective. 
Classical rotating bodies such as galaxies (period=10 15 s), stars (10 6 s), and fluid droplets (Is), and 
quantum rotating nuclei (10 -2O s) may be described in terms of a single 'subgroup GCM(3) (for 
general collective motion in 3 dimensions) of the noncompact symplectic Lie group Sp(3,R). In 
classical physics, the GCM(3) theory is identical to the Riemann ellipsoidal model [1, 2, 3]. 

A Riemann ellipsoid is a uniform density fluid with an ellipsoidal boundary whose velocity field 
is a linear function of the inertial frame Cartesian position coordinates X. The isodensity surfaces 
of elliptical galaxies are very nearly ellipsoidal [4]. Linear velocity fields U L (the superscript L 
indicates a laboratory inertial frame quantity) span the dynamical continuum from rigid rotation, 
IJ L (X) = U L x X, t,o inotational flow, V x () l = 0. Thus, Riemann ellipsoids can model a wide 
class of rotating systems. 

The principal aim of this paper is to present the classical symplectic model with particular em- 
phasis upon its relationship with the Riemann ellipsoidal model [5]. But first the Riemann model 
ami its equivalence to the algebraic CCM(3) theory will be reviewed. To describe a linear velocity 
field, the dynamical group GCM(3) contains the general linear group GL(3,R) as a subgroup. In 
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addition, to characterize the size, deformation, and orientation of an ellipsoid, the GCM(3) Lie 
algebra includes the inertia tensor. 

There are several advantages to adopting the powerful dynamical group method. First, the 
Euler fluid equations of motion for a Riemann ellipsoid can be proven to form a Hamiltonian 
dynamical system [3, 6]. A Riemann ellipsoid phase space is a co-adjoint orbit of GCM(3), and its 
Poisson bracket is inherited from the Lie algebra structure of GCM(3). Moreover, this Hamiltonian 
system is a special Lax pair system [7]. Second, the group method is not restricted to continuum 
fluids. GCM(3) dynamical symmetry applies equally well to discrete systems of particles. Third, 
GCM(3) dynamical symmetry also applies to some quantum rotating bodies. For example, the 
Bohr-Mottelson irrotational surface wave model of collective rotational and vibrational states 
forms an irreducible unitary representation of GCM(3) [8, 9, 10]. Finally, GCM(3) symmetry 
suggests a natural extension to symplectic Sp(3,R) dynamical symmetry [11]. The latter replaces 
the collective kinetic energy of the G0M(3) theory by its exact microscopic expression. 

The hydrodynamic Riemann ellipsoidal model provides a physical interpretation to the abstract 
GCM(3) theory: The length C of the Kelvin circulation vector, a constant of the motion for a 
frictionless, homoentropic fluid flow, is the Casimir invariant for GCM(3) [6]. 

The velocity fields of rigid rotors and irrotational droplets have very different Kelvin circulation 
vectors C. Suppose the rotating system has an ellipsoidal boundary with semi-axes lengths a*. 
The inertial frame Kelvin circulation vector, projected onto the Arth body-fixed axis, is defined as 
the line integral of the velocity field V around the boundary of the ellipse Dk in the i — j principal 
plane for i, j, k cyclic. According to Stoke’s theorem, these line integrals equal the surface integrals 
of the curl of the velocity field, 
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where U denotes the projection of the inertial frame velocity field onto the body-fixed axes, and 
M is the fluid’s masi. 

By definition, the curl of the velocity field of an irrotational droplet is zero, and, hence, 
the Kelvin circulation of an irrotational fluid vanishes, (7 = 0. For a rigid rotor velocity field, 
V x U = 2u;. Because 7r muj is the area of the ellipse D the rigid rotor circulation components 
equal Ck = (2M/b)(ii(iju>k ■ For a general linear velocity field, the curl is a constant vectorfield V x 
(7 = £ - 1- 2u5, where (j is called the uniform vorticity. As the uniform vorticity ranges continuously 
from zero to the negative of twice the angular velocity, the complete Riemann sequence from rigid 
rotation to irrotational flow is traversed. 


2 GCM(3) Dynamical Symmetry 

The symplectic algebra Sp(3,R) consists of the inertia, virial momentum, and kinetic tensors [11]: 

C ij 7 ' Xgj A aj , 

Af = E-V.i/’.i, (2) 

= T. m « 
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where the sums are carried over the particle index a — 1, . . . , A, ra 0 denotes the mass of particle 
o, and X„, P a are the inertial frame vector Cartesian position and momentum of particle a. In 
fluid dynamics, the sums over particles are replaced by integrals over the density distribution, e.g., 
Q\. - f p(X)XiXj(PX . The Poisson brackets close to form the symplectic algebra: 

{^X,} = 

{<&.<&} = 

{$j,Nh} = « ii 

{rJj,74} = o, 

{4 ,7S} = SuTji+SvTf, 

{Q\'j -^ti ) = + GjlNjk- 

The general collective motion GCM(3) subalgebra includes only the inertia tensor Q\ 3 and the 
virial momentum tensor Njj. The* rotational R()T(3) subalgebra is spanned by just the inertia 
tensor and the antisymmetric part oi the viral momentum tensor, viz., the angular momentum 
The Lie algebra CJL(3,R) of the general linear group is generated by the virial 
momentum tensor, and the Lie algebra SO(3) of the rotation group is, generated by the angu- 
lar momentum. The inertia tensor generates a 6 dimensional R 6 abelian Lie algebra. GCM(3) 
and ROT(3) are semidirect sum Lie algebras of the abelian ideal R 6 with GL(3,R) and SO(3), 
respectively. 

In the principal axis frame, the inertia tensor Q is, by definition, diagonal, and its eigenvalues 
are proportional to the squared axis lengths af ot the inertia ellipsoid. 

Although the exact kinetic tensor T tJ is not an element of GCM(3), its linear velocity field value, 
the collective kinetic tensor, is a function ol the algebra generators, [12] t — • Q 1 * N . The 

Kelvin circulation ol a linear velocity field may be expressed in terms of the GCM(3) generators 
asC k = tkiAQ-' ,2 'N*Q l/2 te. 

Time evolution in the classical collective models based upon ROT(3), GCM(3), and Sp(3,R) 
is governed by Hamiltonian dynamics of a special type known as a Lax system. Consider first 
the simple case of ROT(3) for which the dynamics corresponds to Euler rigid body rotation. If 
the inertia ellipsoid is rotating with an angular velocity and L x j = — N JX is the 

angular momentum tensor, then Hamiltonian dynamics is given by 

L = [HL\. (3) 

In terms of vectors, this equation is the familiar law L = -uJ x L that, determines the precession 
of the angular momentum vector in the body-fixed frame. 

A matrix equation of the form X ~ [ F, X ] is called a Lax equation and X — F are referred to 
as a Lax pair [13, 14]. A useful property of any Lax equation is that the trace of any power of X 
is conserved. Let / ;i denote the trace of the pth power of the matrix X . 

l v = -Tr(AT. (4) 

V 
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( 5 ) 


For any Lax equation, it is evident that every l v is a constant of the motion, 

i p = Tr • A') = Tr (A'" -1 • [F, A']) = Tr (X^FX - X P F) = 0. 

In the case of the Euler equation, / 2 = — L * L is the negative of the squared length of the angular 
momentum vector. If p is odd, then I v is zero. If p > 2 is even, then I v is a function of the squared 
length of the angular momentum vector. Thus, there is only one independent invariant among 
the Lax invariants. 

Suppose that -V (4) is a solution to the Lax equation, X = [F, X], corresponding to the initial 
condition X = Xq. If g(t) is a smooth curve of invertible matrices satisfying the matrix differential 
equation g = F • g with the initial condition g = /, then the solution to the Lax equation is just 
the isospectral deformation, 

X(t) = cj(t)-X 0 -g(t)-'. (6) 

This is proven using the identity dg~ x jdt = —g“ l gg“ l * If ft is constant, the matrix differential 
equation g = i\ • g for the Euler equation has the unique solution g(t) = exp(f It) for the initial 
condition #(0) = I. Thus, g(t) is a curve in the rotation group SO(3), and the isospectral 
deformation L(t) = y(/)L 0 g(t)~ l describes explicitly the precession of the angular momentum 
in the body-fixed frame resulting from the rotation g(t) of the intrinsic frame relative to the 
laboratory frame. Because of the choice of initial conditions for </, Lo represents the constant 
angular momentum vector in the inertial laboratory frame. 

To present the time evolution for Riemann ellipsoids as a Lax system, suppose the potential 
energy in the body-fixed frame* V = V'(aj , a 2 , u :J ) is a smooth function of the axes lengths. For 
a star or galaxy, V is the attractive* gravitational self- energy. For a nucleus, V may be approxi- 
mated by the sum of the attractive* surface* energy and the repulsive Coulomb energy. Define the 
Chandrasekhar potential energy tensor W in the rotating frame to be the diagonal matrix, 

dv 

ij ~ &ijdi ^ ( 7 ) 

and, to impose a constraint to constant volume, define the pressure tensor LI = pv to be the 
product of the hydrostatic pressure p times the ellipsoid’s volume v = Anaiaias/S. Hamiltonian 
dynamics for Riemaim ellipsoids is given as follows [7]: 


Theorem. If the inertia ellipsoid is rotating with an angular velocity f] x y = then the 

Riemann ellipsoid Hamiltonian dynamical system is equivalent to the Lax system, X = [F 1 , X], 
where the 6x6 real matrices X and F in the body-fixed frame are given by 


A' = 


N -Q \ F _ / ft 

t — l N )' \ (W + n)-Q~ x 



( 8 ) 


The quadratic Lax invariant, equals the negative of the squared length of the Kelvin circulation 
vector, / - 2 = Ti(N 2 — t • Q) = —(' 2 . The higher order Lax invariants are either zero (odd powers) 
or are functions of the circulation vector s squared length. 

The phase space for a Riemann ellipsoid obeying the Lax equation is a co-adjoint orbit of 

CCM(3): 
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Theorem. Each Riemann ellipsoid orbit is diffeomorphic to some coset space of GCM(3). The 
coset depends upon the value of the circulation C: 


O c 


GCM{i)/SO( 2) = R u x 5 2 , C ^ 0, dim = 14 
GCM{’i)/SO(‘i) = R u , C = 0, dim = 12 


( 9 ) 


The degenerate orbit is diffeomorphic to 12-dimensional Euclidean space. This irrotational 
flow phase space, coordinatized by a 2p , 7r 2;o cv 0 , tt 0 for the quadrupole and monopole degrees of 
freedom, was quantized by A. Bohr. The generic orbits C ± 0 were undiscovered for many years 
because the significant role of Lie groups in this problem was not appreciated by the Copenhagen 
school. The generic orbits are diffeomorphic to the Cartesian product R 12 x S 2 of Euclidean space 
with the two-dimensional sphere. The topology of the sphere forces the circulation to be quantized 
to integer multiples of Ti in a way parallel to the usual angular momentum quantization. Thus, the 
spectrum of the squared length o( the quant um circulation operator is quantized to ( (0 -f* l)ft , 
where C is a nonnegative integer. 


3 Sp(3,R) Dynamical Symmetry 

Classical symplectic Sp(3,R) time evolution ill the rotating frame is given by the Lax equation, 
A' = [F, X], if, in the Lax matrix A\ the linear approximation t to the kinetic energy is replaced 
by its exact expression T. In this way, the restriction to linear velocity fields of the Riemann 
GCM(3) model is removed in the symplectic Sp(3,R) theory. 

The symplectic conservation laws are provided by the Lax invariants I p . The quadratic Casmir 
invariant of the symplectic algebra is the quadratic Lax invariant, C ^ = Tr (N 2 — Q • T). Note 
that for a linear velocity field, the quadratic symplectic invariant simplifies to the negative of the 
squared length of the Kelvin circulation vector. The odd order invariants vanish. The quartic 
symplectic Casimir invariant is the quartic Lax invariant, 

C w = Tr [{NQ - Q l N)(TN - l NT)] - 1/2 Tr [(N 2 - QT) 2 } . (10) 

There is only one more independent Casimir and Lax invariant C (6) = / G ; the higher order invari- 
ants are functionally dependent upon the three independent Casimirs = I p for p = 2,4,6. 

Since the matrices .V and F are elements of the symplectic Lie algebra, the following theorem 
may be proved: 

Theorem. Every solution to the classical symplectic Lax system is given by a isospectral trans- 
formation g{t) 6 Sp(3,R) applied to the initial state 

X{t) = y(t) -X 0 - g(t)-\ (11) 

where X 0 and A" are elements of the symplectic Lie algebra sp(3,R). The group element g(t) is a 
solution to the matrix differential equation g = Fg with the initial condition g = / if and only if 
A" is a solution to the Lax equation with the initial condition X = Aq. 
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( 12 ) 


Consider the co- adjoint orbit of the symplectic group through the point X , 

O x = {g-X-g- 1 | <,eSp(3,R)}. 

The co-adjoint orbit is regarded as a surface in the Euclidean symplectic dual space, sp(3, R)*. A 
manifold that intersect each co-adjoint orbit exactly once is called a “transversal.” A transversal 
T for the symplectic co-adjoint group action is provided by a three-dimensional surface [15, 16] 

T = {.? = ( ® “ 5 ) € sp(3,RY | 5 = diag(s 1 ,s 2 ,s 3 )| . (13) 

Transversal points correspond to elementary systems for which the virial momentum tensor van- 
ishes, N — 0, and the inertia and kinetic tensors are equal and diagonal, Q = T = S. Since the 
inertia and kinetic tensors are positive-definite, the physically relevant transversal consists of only 
those points for which S is positive-definite, .s, > 0. 

An orbit of the transversal point $ E T is diffeomorphic to a coset space of the symplectic 
group modulo the isotropy subgroup. These isotropy subgroups may be proven to be subgroups 
of the unitary group, 

</(3) a | f " ~ v j € Sp( 3, R) | U + iV e 1/(3) ] , (14) 

and, thereby, the coset spaces are given explicitly as follows [15, 16, 5]: 

Theorem. The symplectic phase spaces are diffeomorphic to coset spaces of Sp(3,R): 

.S';/(3, R)j[U{ 1 ) x U( 1) x f/( 1 )], S{ distinct, dim = 18 

Os = | ,S>(3, R)/[U{ 2) x (/(!)], 5, = s 2 ^ s 3 , dim = 16 (15) 

, .S>(3, R)/U{ 3), 5l = s 2 = s 3 , dim = 12 

The degenerate orbit Sp(3,R)/U(3) is diffeomorphic to the complex Siegel half- plane. In future 

work, the dynamical system on the Siegel half-plane will be reported. 
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